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1 Introduction 



A few years ago Chowdhury, Ackerson and Clark [1] reported an interesting light 
scattering experiment on laser induced freezing of a two-dimensional suspension of 
strongly interacting colloidal particles. They showed that the colloidal liquid freezes 
into a two-dimensional crystalline phase with predominantly hexagonal order, when 
it is subject to a one-dimensional potential induced by a standing wave pattern of 
interfering laser beams. The wavevector q Q of the modulating potential was chosen 
to coincide with the ordering wavevector of the colloidal liquid, i.e., the location 
of the first peak of the direct correlation function (DCF) c?'(q). This observation 
motivated experimental studies involving direct microscopic observations [2] and 
Monte-Carlo simulations [3] which confirmed the existence of such a laser induced 
freezing transition. However, the conclusions regarding the nature of the transition 
between the modulated liquid and the crystal have not been definitive. In the orig- 
inal paper where the experimental results were reported, Chowdhury et al.[l] also 
theoretically analyzed the phenomenon of laser induced freezing in terms of a sim- 
ple phenomenological Landau-Alexander-McTague theory and concluded that the 
transition from the modulated liquid to the 2-d (modulated) crystalline phase can 
be made continuous for sufficiently large laser fields. The phenomenological Landau 
theory, based on a polynomial expansion of the free energy in powers of the order 
parameter truncated up to a finite degree, has well-known limitations in the case 
of a first order transition where the order parameter shows a jump. Apart from 
this, the coefficients in a phenomenological Landau type of free-energy are unknown 
and have unknown dependence on experimental parameters. In contrast the density 
functional theory (DFT) [4] pioneered by Ramakrishnan and Yusouff has led to a 
very successful approach for studying the strongly first order liquid-solid transitions 
for various systems including colloidal suspensions [5]. In contrast to the Landau- 
type free-energy functions, the density functional free-energy [4] is a non-polynomial 
function, (important) terms involving arbitrary powers in the order parameters be- 
ing present. Furthermore, the coefficients of the density functional free energy are 
determined from the experimentally measured liquid correlation functions. It is of 
obvious interest to carry out DFT studies of Laser Induced Freezing. Recent density 
functional studies [6,7] have dealt with Laser Induced Freezing transition. 

In reference [7] it has been shown how the freezing (modulated liquid — > crystal) 
transition changes from a first order to a continuous one via a tricritical point with 
the increase in the external potential V e , if the modulation wave-vector is tuned at 
the ordering wavevector of the liquid. The order parameters in this theory [7] are 
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the Fourier components of the molecular field £(r)(= ln(p(r)/p )) with wavevectors 
equal to the reciprocal lattice vectors (RLV) of the periodic structure, into which the 
liquid would have frozen in the absence of V e . Here p(r) is the local density in the 
modulated liquid or the crystal and po is the mean density of the liquid. Symmetry 
considerations show that in the presence of V e , the order parameters corresponding 
to the smallest set of RLVs of the crystal divide into two classes: (1) those corre- 
sponding to the RLVs parallel to the modulation wave- vectors {g/} of V e denoted 
by £/; and (2) those corresponding to the rest of the RLVs in the smallest set {gd} 
denoted by One can choose the wavevectors pertaining to the external potential 
in such a way that an integral combination of vectors in the class {g/} cannot be 
obtained from an odd combination of vectors in the class {gd}. Under this condi- 
tion it is shown [7] that the Landau free energy obtained by expanding the DFT 
free energy about the modulated liquid phase (£/ 7^ 0) in powers of £d has only even 
powers of £d and the transition changes from first order for low values of —j3V e [7] 
(where the quartic coefficient of the Landau free energy is negative) to a continuous 
one for large values of —f3V e [7] (where the quartic coefficient is positive). Note that 
this symmetry argument holds good under the inclusion of arbitrarily large RLV 
and liquid direct correlation functions of any order [7]. 

The analysis of reference [7], however, is of a mean field nature. We clarify below in 
what precise sense the theory presented in reference [7] is a mean field theory. The 
density p(v), or equivalently the molecular field £(r), contains the freezing modes, 
i.e., the modes with wavevectors equal to the reciprocal lattice vectors (RLV) of 
the periodic structure (into which the liquid would have frozen in the absence of 
V e ), and other modes which we call the fluctuating modes. So we can write that 
£ ( r ) = S g £g ex P(^g ' r ) + £q exp(zq • r) where the first summation runs over the 
freezing modes and the second summation runs over the fluctuating modes. In the 
present problem, for instance, £/ and £d are the freezing modes corresponding to the 
smallest RLVs parallel to the wave-vectors of V e and the rest of the RLVs in the 
smallest set respectively. In reference [7] only the freezing modes are retained and 
all the fluctuating modes are completely neglected. Thus the theoretical analysis 
presented in reference [7] is a mean field one. In general, the fluctuating modes can 
have an important effect on the freezing transition. This motivates us to extend the 
mean field treatment of reference [7] by including the fluctuating modes in addition 
to the freezing modes which are already present in the mean field study of reference 
[7]. In this context we recall the work of Brazovskii [9-11]. 

He considered the transition from an isotropic liquid to a cholesteric liquid crystal. 
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Here the mean field Landau free energy in terms of the order parameter contains 
only even powers and T^ A \ the coefficient of the quartic term is positive; hence the 
transition is predicted to be continuous. He demonstrated that fluctuations trans- 
form a continuous transition to a first order one. The large effect of fluctuations 
in Brazovskii's theory is due to a considerable softening of the order parameter 
modes at a nonzero value qo of the wavevector. As a result, the effect of fluctuations 
comes from a large surface area (~ ^irq"^) in the momentum space. After integrat- 
ing out these fluctuation modes, the coefficient of the quartic term in the effective 
Landau free energy becomes negative where the coefficient of the quadratic term 
changes sign, indicating that the continuous freezing transition is preempted by a 
first-order one. So an obvious point to investigate is to what extent the fluctuating 
wavevectors close to the ordering wavevector q Q of the liquid are important in the 
presence of an external modulation — (3V e tuned at go [1-3,7]. More specifically we 
ask if the quartic coefficient in the effective Landau expansion can become nega- 
tive (where the quadratic coefficient changes sign) in the context of Laser Induced 
Freezing transition in the large — f3V e limit, where the transition is predicted to be 
continuous in reference [7]. We consider a two-dimensional colloidal system subject 
to a one- dimensional external potential for the sake of simplicity. Furthermore, it is 
well-known that the effect of fluctuations is more prominent in two dimensions than 
in three dimensions. 

Our starting point is the free-energy cost of creating a density inhomogeneity over 
a uniform liquid (in the incompressible limit), in terms of £ q [7,8]: 
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(i) 



where the summation over q includes the set ({gd}) and the fluctuating modes. Here 

(2) 

c q is the direct correlation function of the liquid, V is the volume of the system and 
rif is the number of vectors in the class ({gf}). We generate a Landau expansion 
from equation (1) about the modulated liquid phase (£f ^ 0) in powers of £ q (where 
the transition is continuous according to the mean field theory) and then integrate 
out the fluctuating modes to form an effective Landau expansion in powers of £d- 
This free energy has been used to obtain the phase diagram. Note that the odd 
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invariants of £d are forbidden in the effective Landau free-energy due to symme- 
try restrictions as in the mean field case [7]. In this work we follow Brazovskii's 
approach to form the effective Landau free-energy. In particular, Brazovskii shows 
that the most dominant corrections to the two-point and four-point correlation func- 
tions come from diagrams with a single loop. As a consequence, it is sufficient to 
evaluate self consistently the self-energy corrections to one-loop order. Similarly, it 
is shown [9] that the corrections to the four-point vertex functions are dominated 
by a restricted class of one-loop diagrams and their ladders. Here we have a liquid 
of strongly interacting colloidal particles characterized by a static structure factor 
S(q), with a sharp maximum at q ~ q . Equivalently, this implies that [^(g)] -1 
has a prominent minimum at q ~ qo, a feature present in Brazovskii's theory as 
well. However, in this situation there are a few additional complications which are 
absent in Brazovskii's problem. Here we consider a phase transition induced by an 
external potential. The coupling of the external potential to certain density modes 
of the isotropic liquid leads to a nontrivial structure of the correlation matrix in the 
modulated liquid phase which plays a central role in the computation of the renor- 
malized four-point vertex functions. Furthermore, the momentum conservation laws 
appearing in the present work involve the additional wavevectors pertaining to the 
external modulation. Also, the bare vertex functions appearing in our free energy 
have a complicated dependence on the external potential. For the sake of simplic- 
ity we confine ourselves to asymptotically large values of the external potential. 
In this limit, we neglect fluctuations in the modes £f (pertaining to the external 
potential) since a zero temperature normal mode analysis shows that the external 
potential creates a gap (A ~ \Ve\) in the phonon spectrum along the direction of 
externally induced ordering. Consequently fluctuations in these order parameter 
modes are energetically unfavorable for \/3V e \ ^> 1. We therefore treat the modes £f 
non-perturbatively [7]. We calculate the correlation matrix in the modulated liquid 
phase by considering one-loop Hartree corrections to the self-energy (Sect. 3). It is 
easy to evaluate the coefficient (T 2 ) of the quadratic term in this effective Landau 
free energy from a knowledge of the correlation matrix in the modulated liquid phase 
(Sect. 3). We determine the point in the phase diagram where T 2 goes to zero. We 
then use the correlation matrix to evaluate the renormalized four-point vertex func- 
tion and consequently T4, the coefficient of the quartic term in the effective Landau 
free energy at that point (Sect. 4). We find that T4 remains positive in our regime 
of interest indicating that the freezing transition remains continuous for large values 
of V e even in the presence of fluctuations. We also show, as we will elaborate later, 
that fluctuations tend to stabilize the liquid phase relative to the solid phase in the 
limit of large V e . 
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The paper is organized as follows. In Sections 2, 3 and 4, we give a detailed account 
of our calculation leading to the determination of the nature of the laser induced 
freezing transition and in Section 5 we present our numerical results. We conclude 
the paper in Section 6 with a few remarks. 



2 Free Energy 

Expanding the free-energy in equation (1), in a power series and retaining terms to 
quartic order in £ q (q 7^ gf or zero), we get: 



-^qi ,Q2 ,q3 £qi £q2 £qa 

q qi,q2 qi,q2,q3 

_ 9/L ^ / Qqi,q2,q3,q4^qi^q2^q3^q4 (2) 
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qi ' q2 ' q3 /d^expfcV.e^) 



■*3 

Qqi,q2,q3,q4 = < 5q 1 ,q 2 ,q3,q4 + <2qi .q2.q3.q4 

(3) 



with 



/ d d r[£f ^ . e( is f 3)r )]e^^' i = 1 ' 4q ^ 



(i) 

qi,q2,q3,q4 r r . _ /. (j) 



fd'Vexp&^.e^'^ 
'5qi' ) ,q2,q3,q4 = — 3A qi q2 A q3 q4 . (4) 

One can easily verify that each of these coefficients is nonzero if and only if the sum 
of the q-vectors appearing in the coefficient is Gf, which is an arbitrary integral 
combination of vectors {gf}. This is the momentum conservation condition for the 
present problem. For instance, A mm in equation (3) will be nonzero if qi + q2 = Gf. 
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The cubic term in equation (3) will have nonzero value if Yli=x 3 q« = Gf. Similarly 
< 5qi ) q2 , q3 , q4 in equation (4) is nonzero if ^ i=M q* = G f . 

The mean field calculation shows that in the absence of V e a 2-dimensional liquid 
freezes into a triangular lattice [1, 12] with its 6 smallest reciprocal lattice vec- 
tors (RLV) {g } = (0,±l)g 0> (±l/2,±>/3/2)go- Here *>> as is known from DFT, 
corresponds to the first peak of DCF. We choose the wavevectors of V e to be 
{gf } = (0, ±1)50 m this set. This choice simplifies the geometry of the problem 
greatly. With this choice, the vector Gf can be written as mgf where m is an arbi- 
trary integer and gf = (0, l)qo- Since there is a one-dimensional ordering created by 
V e , there is a periodicity in the x-direction. This symmetry of the problem enables 
us to write a general wavevector p as a sum of q, a vector in the first Brillouin 
Zone (BZ) and an arbitrary number of gf. Here the first BZ is the strip between 
q y = —0.5 to q y = 0.5. Now we introduce some simplifications in our notation. Let 
us consider, for instance, ^4 qi , q2 defined in equation (3). Note that the momentum 
conservation demands that qi + q2 = mgf. So if we choose a vector q in the first Bril- 
louin Zone such that qi = q + igf, then the momentum conservation demands that 
q2 = -q + jgf, so that i + j = m. So we define A qiiq2 = ^(q+; gf ),(-q+jg f ) = A ij (q L ), 
where the momentum conservation law £\ q^ = mgf, is automatically satisfied. One 
can write similar simplified notations for the cubic and the quartic vertices. 



In the large — f3V e limit, the integrals appearing in the various coefficients in equation 
(3) can be evaluated by means of an asymptotic expansion. Let us consider the 
vertex in equation (3). For the 2-d hep lattice and the given choice of gf, we have 
<pf = J2i ex P(^gf*' ) r ) — 2cos(^x). In the large — (3V e limit one can expand <pf about 
its maximum and the resulting Gaussian integral can be computed exactly. In this 
way the following expressions can be derived: 



Ai 



T 

J- 1 



o (1) 
o {2) 



exp 



exp 



exp 



3 exp 



Hi + 3) 
647T 2 £ f 

3(i + j + k) 2 ~ 

64vr 2 «ef 

3(i+j + k + lf 

647T 2 Cf 

Hi + jf 
64vr 2 ^ f 



exp 



3(k + l)' 

647T 2 ^f 



(5) 



In the large —/3V e limit the vertices do not depend on the momentum labels. Note 
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that — j3V e enters the calculation through £ f . In equation (1) we set £ q = for all 
q except q = gf and make an asymptotic expansion as above. The self-consistency 
condition at the minimum of the resulting free energy with respect to £f yields the 
following expression for £f: 

& = -/3K + Poc (2) (g ) + l/(W e ) 

to 0(1/ pV e ) in this limit. 



3 Correlation Matrix 

The first step is to calculate the bare correlation matrix for the modulated liquid 
phase. Since the effect of the external potential is to scatter a mode with momentum 
q into that with 




(a) 




Fig. 1. Lowest order correction to the self-energy due to (a) the 4-point vertices 
and (b) the 3-point vertices. 

momentum q + ragf , where m is any integer, the quadratic part of the free-energy 
equation (2) has off-diagonal couplings, the strength of which is given by in 
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equation (5). The elements of this matrix are: 

(diagonal) 

(off — diagonal) (6) 

The bare correlation matrix Q^' for the modulated phase is obtained in terms of 
S- 2 \q) and — (3V e (via £f) by numerically inverting the matrix in equation (6). The 
main difficulty is that the correlation matrix is of infinite order. 

At this point we introduce some approximations that enable us to deal with a 
sufficiently simple and finite dimensional correlation matrix. We estimate the con- 
tributions of the corrections coming from the cubic and the quartic vertices. To the 
lowest order the correction due to the quartic vertex is given by: 

E = E / d P£- (P) • ( ? ) 

ij ^ ' m,n 

This correction does not depend on the external momentum label q (see Fig. la). 
Similarly, the lowest order contribution due to the cubic vertex is: 

F 1 f 

E(q) = E T ihm 2 T jhmi J dpQlX (p)^ m2 (q - P ). (8) 

ij h ,h ,mi ,m2 

This correction, however, depends on the external momentum q (see Fig. lb). The 
integrations in both cases are over the first BZ. We determine the more dominant 
of these two self-energy corrections. Near go we take a simpler form for the liquid 
direct correlation function: 

[l/p c (2) (|q|)]-l = c + %-go) 2 , (9) 

where cq and d parametrize the liquid direct correlation function. 
Note that: 

p c (2) (go) = l/(l + c ), 

i.e., Co pertains to the first peak height of the liquid direct correlation function and 
d corresponds to the width]. We compute Ylfj (Eq. (7)) and (Eq. (8)), using the 



p c( 2 H|q + «gf|) 

( 3(z+jH 
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bare correlation matrix (see Eq. (6)). We go to the diagonal basis and compare 
the eigenvalues of Y2h an d Y2jj f° r var i° us <!• It turns out that the most dominant 
eigenvalue of ^^-(q) gives smaller corrections than that of Ylfj f° r an IqI ly m § m 
the region near go- This means that the predominant correction comes from Ylfj 
for |q| ~ g . We also observe that the most dominant eigenvalue does not change 
significantly with increased dimensionality of Y^fj beyond 3(i = ±1 and j = ±1). 
To summarize, we confine our calculations to the region |q| ~ q (i.e., we confine 
to fluctuations in the low energy modes of the order parameter spectrum), keeping 
only Y^fj m the self-energy correction and restrict to the block E (—1, 1) of Y2ij- 

Notice that the off-diagonal elements of Qq 1 depends only on the magnitude of 
and not on % and j separately. This means that there are only two independent 
off-diagonal elements. Hence within our approximation the 3x3 bare correlation 
matrix is as displayed below: 

c + d(\q - flff| - g ) 2 a b 

a c + d(q - g ) 2 ao 

b ao c+d(\q + g { \ - g ) 2 _ 

Here ao = — exp[— 3/647r 2 £f] is the off-diagonal coupling with (i + j) = 1. Similarly, 
bo = — exp[— (12/647r 2 £f)] is the off-diagonal coupling with (i + j) = 2. We refer to 
c , a , b and d as the bare parameters. 



Go l = 



Dyson's equation for the corrected correlation matrix is: 



(10) 



where 



E(i) 



H 

E 



2(2vr) 



dp^ m n(p) 



'ID 



Since the correction Y^^ does not depend on q, the q-dependent term in the bare and 
the corrected correlation matrix are the same (i.e., d is unaffected by the correction). 
The bare parameter Co gets renormalized to c as a result of the correction and so do 
the off-diagonal elements. Hence, Q~ l is a 3 x 3 matrix: 

c + d(\q-gi\ -g ) 2 a b 



a 
b 



c + d(q- g ) 



c+d(\q + g { \ -q o y 
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where now a and b are renormalized off-diagonal couplings. Hence the self-consistency 
equations derived from equation (10) involve three nontrivial parameters c, a and b. 
It is clear that in order to frame these equations it is sufficient to consider diagrams 
pertaining to momentum transfers (i.e. (i+j)) of magnitude 0, 1 and 2 only. These 
contributions are denoted by 53 (oo) (^ ee Fig. 2), 53(-io) anc ^ 53(-i i) respectively. 
The self-consistency equations 



^(0,0) 



c = c - 


E. 




(0,0) 


a = do — 


£. 




(-1,0) 


b = bo- 


E. 












+ .O. 


it I"' <T 


7 «.<r 



q- (3).<f 3 < 4 > -er 



q- (5) --5 

Fig. 2. - The diagrams which contribute to 53(oov 



(12) 
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where E(o,o)> £)(-i,o) and E(-i,i) are S iven b y : 

E = / d P[^(^0,o(p) + 01,! (P) + 0-l,-l(p)) + <S£l, (p) + ^-1,1(P)] 
(0,0) 

E = / d p[^(^o.o(p) + ^i,i(p) + 0-i.-i(p)) + m, (p) + we-i,i(p)] 

(-1,0) 7 

]T = /"dp[M(0o,o(p) + 0i I i(p)+0-i I -i(p)+^ 1)O (p)+P^i 1 i(p)].(13) 

(-H) 

Here, 

i? = 1/(2tt) 2 

5 = 4exp(-3/647r 2 £ f )/(27r) 2 

T = 2exp(-12/647r 2 ^ f )/(27r) 2 

U = exp(-3/647r 2 ^ f )/(27r) 2 

1/ = (-2 exp(-12/647r 2 £ f ) + 12 exp(-6/64vr 2 £ f ) + 3 exp(-l/647r 2 £ f ) - 2)/(2tt) 2 

W = (-exp(-27/647r 2 ^ f ) + 6 exp(-15/647r 2 £ f ) - exp(-3/647r 2 £ f ))/2(27r) 2 

M = exp(-12/647r 2 £ f )/(27r) 2 

N = (-2 exp(-3/647r 2 ^ f ) + 12 exp(-15/647r 2 £ f ) - 2 exp(-27/647r 2 ^ f ))/2(2vr) 2 

P = (-4exp(-24/64vr 2 £ f ) -exp(-48/64vr 2 £ f ) - 1)/2(2tt) 2 (14) 

Let us consider diagram (4) in Figure 2 for the purpose of illustration. The sum of 
the momenta at the vertex is — gf. Hence the contribution of is exp(— 3/647r 2 £f). 
To calculate the contribution from Q^ 2 \ we note that the sum of the momenta of 
the two free legs is zero and that of the internal lines is — gf. So the contribution is 
— 3 exp(— 3/647r 2 £f). Adding the two contributions and multiplying the sum by the 
symmetry factor and ^2 coming from the phase volume, one gets the expression for 
S in 5^(0 0)' Similar consideration yields the other numerical coefficients. 

We solve the coupled equations in equation (12) numerically for c, a and b for a given 
Co and — (3V e (i.e., given bare parameters cq, ao and bo)- From the self-consistent 
values of c, a and b, one has to construct the coefficient of the term (£d) 2 in the 
effective Landau free-energy. To this end we enumerate the freezing modes as follows: 
gf = (V3/2)(l,l/V3), gf = (V3/2)(l,-lA/3), gf = (V3/2)(-l,-l/V3), 
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as 



pi 



(a) 



>otcxx + >oocx 



(Q 



Fig. 3. - (a) The choice of the external legs for a typical 4-point vertex that appears 
in the quartic term in the effective Landau free energy, (b) A one-loop correction 
to a 4-point vertex with external legs as in (a). The labels on the lines indicate the 
number of gf to be added to the momentum concerned, for instance label l\ indicates 
that the momentum is -Migf. Similarly the labels n\ and n 2 on the internal line 
indicates the ni^th element of the Q matrix, (c) The ladder corresponding to (b). 

gd = (v^3/2)(— 1, l/\/3)- If we assume that all the freezing modes are degenerate, 
the coefficient of the quadratic term is given by: 

T2 = GoM^) + 2^i(g d 1} ) + + Go,M ] ) + Mo,M ] ) + ^(gf )■ 

(15) 

We find the point in c — (—/3V e ) plane for which T 2 = O. 



4 Renormalization of the 4-Point Vertices 

In Section 3 we locate the position at which T 2 goes through zero in the c — (—j3V e ) 
parameter space. Our aim here is to investigate the sign of the coefficient T 4 of 
the quartic term in the effective free energy at this point. The first step towards 
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investigating this question is to evaluate the renormalization of the four-point ver- 
tices appearing in the free energy expansion given in Section 2. Here we outline the 
calculation of the corrections to the four-point vertex functions. 

The terms quartic in the freezing mode £d, correspond to diagrams with four ex- 
ternal legs labelled by the wavevectors belonging to the class {gd} (See Fig. 3a). 
Given the pair of wavevectors g d ^ and gjj = — g^, one can generate the remain- 
ing two vectors in the class {gd} by adding to them a suitable integral multiple 
mgf with m taking integral values between —1 and +1. We restrict ourselves to 
the lowest order {i.e. one-loop diagrams and their ladders) fluctuation corrections 
coming from intermediate interactions between various modes. We use the renor- 
malized correlation functions obtained via Dyson's equation (Sect. 3) to evaluate 
the one-loop corrections to the four-point vertex functions. We notice that for the 
four-point vertices shown in Figure 3a there are no one-loop corrections coming from 
interactions mediated by cubic vertices. This is due to the absence of such cubic 
vertices in the effective free energy which follows from the fact that we choose the 
wavevectors pertaining to the external potential in such a way that an integral com- 
bination of vectors in the class {gf} cannot be obtained from an odd combination 
of vectors in the classjgd}. Thus, corrections to the bare four-point vertex func- 
tions come entirely from intermediate scattering processes involving the four-point 
vertices. We emphasize that the four-point vertices have two parts, and 
(see Eq. (4)). Clearly Q^ 2 ' corresponds to a more stringent momentum conservation 
condition compared to Consider for instance a general one-loop diagram of 

the form shown in Figure 3b. It is easy to see by inspection that corrections to 
come from such a loop only if q = . In contrast, one-loop corrections to 
come from the entire range of values of q. 

We define a function Xm 1 m 2 m / 1 m' 2 as follows: 



+i 




+ n (2) 



m i m 2 n i n 2 ^ n i n 2 rr4 



m=— 1 7i2=— l 



where 
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and 

tt(2) 1 r , (l)\ n ( (1)\ 

Summing over the ladders of one- loop diagrams (Fig. 3c), we obtain the correction 
to the vertex Q mxm im' x m' 2 - 

ni,n 2 

where n — -+- TT (2 ^ Thus 

wneie Vmim 2 nin 2 — ^mim2hki LL l 1 l' 1 kik' 1 ^l' 1 k' 1 n' 1 n' 2 ^ ^mim2hki LL hl' 1 kik' 1 ^l' 1 k' 1 n' 1 n 2 - LLmb 

the corrected four-point vertex function corresponding to external legs labelled by 
mi, m 2 , m'j and m' 2 is given by 

Qm 1 m,2m' 1 m 2 ~ Qmim2m' 1 m' 2 + (^^)mim2m' 1 mj 

with Qm 1 m2m' 1 m' 2 the bare four-point vertex. 

The relevant quantity which dictates the nature of the phase transition is, of course, 
T 4 which is the coefficient of the term quartic in £d in the effective free energy. It is 
given by 



^4 — / , Qm-imim',™!, (16) 

where the summation is over those values of m\ 1 m2 1 m' x and m' 2 which correspond 
to the vertices with external legs labelled by the vectors belonging to the class {gd}- 



5 Numerical Results 

In Figure 4, we show the renormalized c and the bare Cq(= l/poc^ 2 \qo) — 1) for 
different values of — (3V e at which T 2 = (Eq. (15)). We notice that the values 
of c are always larger than cq. This indicates that fluctuations tend to reduce the 
strength of the correlation in the system. The enhanced effect of fluctuations with 
increasing — (3V e may probably be attributed to the partial (1 — d) ordering caused by 
the external potential. We find that T 4 (evaluated from Eq. (16)) at the point T 2 = 
is positive, which is the signature of a continuous transition. So, the transition is 
continuous as found in the mean field calculations as well. 
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Fig. 4. - The values of c (dashed line) and Cq (solid line) at T 2 = for different 

-pv e . 
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Fig. 5. - The phase diagram, i.e., T 2 = line in the poC^(q ) — (—{3V e ) plane 
obtained from the calculations including fluctuation effects (solid line) and from the 
mean field calculations (dashed line). 
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The phase diagram in the poc^ (qo) — (—j3V e ) plane is shown in Figure 5. Notice that 
as one goes to higher —0V e , in comparison to the mean field theory the fluctuations 
enhance the stability of the liquid phase relative to the crystal phase but the transi- 
tion line eventually saturates as in the mean field theory However, in the mean field 
phase diagram the critical line asymptotes to poc^(qo) = 0.5 for large — (3V e from 
above, while that in the present calculation asymptotes to /0o c (?o) — 0.509 from 
below, a feature found in recent simulations as well [13]. Thus there is a difference 
in the curvature of the critical line in the poc^^o) — ( — (3V e ) plane in the two cases, 
which represents a qualitative difference between the present phase diagram and the 
mean field phase diagram. 



6 Conclusion 

In conclusion we have presented a theory of laser induced freezing which accounts 
approximately for the effect of fluctuations on the freezing transition. The freezing 
transition in our calculation remains continuous even in the presence of fluctuations 
for large values of the external modulating potential. This implies that one can 
perform light scattering experiments to look for critical opalescence indicating the 
presence of divergent static correlations of density fluctuations. 

The main qualitative result of this study is contained in Figure 5. As mentioned in 
Section 5, the curvature of the critical line in the poc^ 2 \qo) — (— (3V e ) plane changes 
sign when fluctuations are taken into account. Thus, fluctuations tend to enhance 
the stability of the liquid phase relative to the crystal phase. It would be interesting 
to find out if this feature of the freezing transition is seen in real experiments. 

Due to the inherent complexity of the problem, we have confined ourselves to self- 
consistent one-loop corrections to the self-energy and incorporated only the one-loop 
diagrams and their ladders into the renormalization of the four-point vertex. Fur- 
thermore, we have confined ourselves to the dominant corrections to the vertices 
coming from the low energy region (i.e. |q| ~ q ) of the order parameter spectrum. 
At this level of approximation there are no significant corrections to the self-energy 
from the interactions mediated by the cubic vertices. However, for small values of 
q, these corrections can be quite significant (compared to the one-loop correction 
mediated by the quartic vertices included in this paper). In future it would be cer- 
tainly worthwhile to carry out a more systematic analysis of the role of fluctuations 
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in the laser induced freezing transition. 
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